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FOREWORD 

This report was prepared by The University of Michigan on 
National Aeronautics and Space Administration Task Order Contract 
NASr-5U( lo) , Navigation Satellite Studies. The work was adminis- 
tered under the direction of Mr. Eugene Ehrlich, Chief, Navigation 
and Traffic Control, Space Applications Program Office, National 
Aeronautics and Space Administration. 

The work reported covers the quarter 1 April through 30 
June I 967 * The work was performed by the Navigation and Control 
Systems Group of the Radio Science Laboratory, one of the Willow 
Run Laboratories of The University of Michigan’s Institute of 
Science and Technology. 

The originator’s report number is 7657-12-P. 
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SUMMARY 


This progress report covers the period from 1 April I967 
through 30 June I967. A primary objective of this NASA contract is 
to examine the feasibility of various means of fixing aircraft and 
ship positions through the use of satellites. A major effort on 
this contract has been to generate a general position-error analy- 
sis suitable for examining satellite navigation systems through 
the use of digital computer models. 

This report describes efforts during the current quarter 
regarding aircraft positional fixes using the VHF channels on 
the currently operational satellite ATS -1 and the soon to be 
launched satellite ATS-C, The digital computer model is being 
modified to include VHF propagation effects 9 and the model will 
be used to predict expected aircraft position fixes. In addition, 
equipment characteristics for use aboard aircraft are being 
delineated . 

No effort was expended on the interferometer measurement schemes 
in this quarter. 
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INTRODUCTION 

This progress report ^ covering the period from 1 April 1967 
through 30 June 1967, presents the results obtained from the work 
conducted under NASA Contract NASI 4 - 5 MIO). A primary objective of 
this contract on satellite navigation studies is to examine the feas- 
ability of various means of fixing aircraft and ship positions 
through the use of satellites. To this end, the previous year r s work 
generated a general posit ion-error analysis suitable for examining hy- 
perbolic and range position-fixing systems. Several systems were ex- 
amined . 

During the current quarter, the effort was shifted to the task of 
adapting the model to determine how well aircraft position fixes could 
be obtained when ranges are measured to two geostationary satellites. 
In particular y the Applications Technology Satellites (ATS), ATS-1 and 
ATS-C, both equipped with VHF transponders } were considered. Aircraft 
positions are to be determined from the satellite positions y range 
measurements to the satellites ^ and the aircraft altitude. 

An experiment conducted with these satellites can take many forms 
in terms of equipment on the ground and in the aircraft. To minimize 
the time required to set up an experiment and to minimize cost, it 
would be desirable to use existing facilities and equipment to the 
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greatest practical extent. 


Any practical experiment will consist in either of two forms: 
one in which range measurements to the satellites are made aboard the 
aircraft and one in which range measurements are made at a ground 
station (or stations); the ranges measured are those from the aircraft 
to each satellite and from the ground station to each satellite. 

Position-fix errors arise from imprecise knowledge of satellite 
positions , aircraft-altitude errors , and range-measurement errors. 

But for geostationary satellites , their range positions can be well 
predicted and any error due to satellite motion can be minimized. Also, 
previous calculations on the computer model have shown the effect of 
aircraft-altitude error on position fixes is small , that is to say, 
and altimeter error is well within the tolerable bounds allowable in 
determining position fix. The major effect then in calculating air- 
craft-position errors is caused by the range measurements between the 
satellites and the aircraft. One cannot, however, disregard the other 
errors including those mentioned above and any timing errors arising 
in satellite or aircraft transponders. 

Several error sources exist in measuring range. Precise time 
measurements form the basis for accurate range measurements. Vari- 
ations in transponder delays, whether on board satellites or aircraft 
and any deviations from straight lines in propagation paths will cause 
errors in range measurement. Propagation errors will be the most 
significant ones at VHF. Refraction of the radio waves, both in the 

ionosphere and atmosphere, will cause a propagation path to a satel- 
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lite to lengthen , with the ionosphere causing the greater increase in 
path length. 

Refraction effects are being modeled simply. If sufficient sta- 
tistical data were available on the refractive index of the ionosphere 
on a diurnal , seasonal , and geographical basis y it would be possible to 
obtain a distribution of range errors for given times and for given pene- 
tration angles of the ionosphere by the radio paths. 

A few calculations with the model have yielded values which indi- 
cate path lengthening of about 1000m. The lengthening can be attributed 
almost completely to the ionosphere , especially for the range measured 
to an aircraft flying at high altitudes. 

Note that a model such as this can be made to reduce the magnitude 
of the range -measurement error. But in determining aircraft-position 
fixes , the same model cannot be used to determine and then to correct 
for errors in range measurements arising from propagation conditions. 
Consequent ly one has the choice of specifying a set of conditions as 
those actually existing and then calculating corrections with respect to 
that set of conditions, or one can examine and use data on range-measure- 
ment errors obtained from an experiment. The former course is being 
pursued at present; efforts will be made to obtain the requisite infor- 
mation from the ranging experiments being made with the ATS-1. Also, 
suggestions will be made as to how the effects of propagation conditions 
might be minimized during the course of the position-fixing experiment. 

The following section describes the propagation model as it pre- 
sently is constituted. 
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TROPOSPHERIC AND IONOSPHERIC EFFECTS ON 
RADIO RANGE MEASUREMENTS 

Errors in the measurement of distance made with radio equipment 
are caused in part by the refractive nature of the troposphere and ion- 
osphere. The velocity of propagation of an electromagnetic wave is a 
function of the refractive index over the propagation path, and the path 
itself is distorted from a straight line by refractive bending of the 
wave front. Errors thus incurred in measuring range by assuming straight- 
line propagation at a velocity of light in a vacuum are called radio 
range errors. These errors are of considerable importance in determin- 
ing the accuracy of navigational systems which employ measurements of 
range, range differences, range rate, or any combination of these quali- 
ties . 

The actual value of the range error occurring over a particular 
propagation path is a function of the structure of the index of refrac- 
tion at the particular time that the range measurement is made. The 
refractive-index structure is always to a great extent an unknown quan- 
tity; however, it has a spatial and a time variation. A refractive 
model that is usually employed is one in which the actual refractive 
index is approximated with a structure which is a function only of 
height above the surface of a spherical earth. With such a model, ray 
tracings can be performed and the average effects of the refractive-in- 
dex range errors can be analyzed. 

In the section that follows, an approximate method will be described 
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to calculate the range error due to tropospheric and ionospheric re- 
fraction. 

2.1. COMPUTATIONAL METHOD 

The method presented here is based on that of Weisbrod and Ander- 
son.* This method is particularly simply and can be applied to both 
tropospheric and ionospheric bending. The only assumptions are that 
the refractive-index gradient is only in the vertical plane , that the 
index of refraction can be approximated by a number of linear segments, 
and that the thickness of these steps is small compared to the earth’s 
radius. In this method, there is no limitation on the shape of the 
profile or the angle of elevation. These assumptions are justifiable 
for the problem under consideration. 

2.2. RANGE ERROR IN THE TROPOSPHERE 

The signal retardation dx caused by a layer of thickness dp 
(fig. l) is given by 

d T — (~ — l)csc p dp 
or 

d T = (n-1 ) ~ - £ dp 
c 

where c and v are the signal velocities in free space and the medium, 

Q 

respectively, and n = — . The range error is given by 


Ar 


,p. 


jk 


c dx 


P J 



(n-1 )csc 


p dp 


( 1 ) 


S. Weisbrod and L. J. Anderson, "Simple Methods for Computing 
Tropospheric and Ionospheric Refractive Effects on Radio Waves," Proc . 
I.R.E., Vol . 1+7, October 1959, PP- 1770-1777- 
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By defining 


( 2 ) 


N = (n-1 ) x 10 6 
one can rewrite the above equation as 


Ar 


jk 



N y IQ 
sin p 


( 3 ) 


This integral can be approximated by the expression 


Ar jk - (\ + N j> < p k - P j> X 10 
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sin ft + sin p . (4) 

k J 

Comparison of the exact solution with (4) indicates an error of a few 
hundredths of one percent for the case of the tropospheric propagation. 
In the ionosphere the errors are larger, but sufficiently small to 
justify the use of the approximate expression. The total error is given 
by 


Ar= E 

k=o 


("k + l + V < P k + l - P k> X 10 

(sin P k+1 + sl " p k ) 


If the ray departs from the earth 1 s surface with an elevation 

angle a , Snell’s law for spherical stratification is 
o ’ 

n a cos a = np cos B = constant 
o o 

where n = surface index of refraction 
o 

a = earth’s radius 
p = a + h 

n = index of refraction at specified height 
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From Snell’s law. 


sin p 


n a 
o 


np 



1/2 


2.3* RANGE error in the ionosphere 


(5) 


It can be shown that the index of refraction in the ionosphere, 
excluding the effect of the magnetic field and collisions, is given by 

Ne e2 

n = 1 - 

2 

me ocd 


or 


n 


1-b 


Ne 


(X>2 


where b = e2 
me o 

Ne = electron density per cubic meter 

m = electron charge 

eo = permittivity of free space 

0 ) = angular frequency of the incident signal 

To compute the ionospheric error, one can use equation k, pro- 
vided that the absolute value of N is used. Thus, 


n-1 


Ar 


■L 


k=o 


( N 


k+1 


+ N, 


) (P k+1 - p k ) X 10 


-6 


( sin P k+1 + sin pj 


( 6 ) 


By combining the two errors, one obtains the total range error 
for oneway passage through the troposphere and the ionosphere. 
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